Electrothermal flow in Dielectrophoresis of Single-Walled Carbon Nanotubes 
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We theoretically investigate the impact of the electrothermal flow on the dielectrophoretic separa- 
tion of single- walled carbon nanotubes (SWNT). The electrothermal flow is observed to control the 
motions of semiconducting SWNTs in a sizeable domain near the electrodes under typical experi- 
mental conditions, therefore helping the dielectrophoretic force to attract semiconducting SWNTs 
in a broader range. Moreover, with the increase of the surfactant concentration, the electrothermal 
flow effect is enhanced, and with the change of frequency, the pattern of the electrothermal flow 
changes. It is shown that under some typical experimental conditions of dielectrophoresis separation 
of SWNTs, the electrothermal flow is a dominating factor in determining the motion of SWNTs. 

PACS numbers: 47.65.-d,82.45.-h,73.63.Fg 



I. INTRODUCTION 

Single-walled carbon nanotubes are key materials in 
nanotechnology as potential candidates for diverse apph- 
cations owning to their extraordinary mechanical, ther- 
mal, optical and electrical propertied. On explor- 
ing the utility of the electrical properties, one of the 
current critical challenges is the separation of metallic 
(m-SWNTs) and semiconducting SWNTs (s-SWNTs). 
Among various post-synthesis separation methods de- 
vised and applied-^ 3 4 a 6^ dielectrophoresis separation has 
been demonstrated to be possible with high selectivity 
and simplicity^. Furthermore, DEP also allows us to 
deposit SWNTs to selected sites to construct electric 
device^. While the complexity in the optical measure- 
ments gives rise to difficulties in interpretation of the ob- 
tained spectrsP, deeper understanding in the transport 
dynamics of SWNTs under DEP operation is of a great 
importance. It also allows us to discuss the efficiency 
and possibility for optimization of system designs and 
working conditions. Here, numerical simulations should 
serve as a powerful tool since the dynamical observation 
of SWNTs transported in the suspension is extremely 
challenging in experiments. 

On modeling a SWNT as a prolate ellipsoid, the DEP 
force on the SWNT in an inhomogeneous AC electric field 
E is calculated 



Fdep - ^-^aV\n\ 



(1) 



where a = Re[ J " ], e* = e,n,p-i^ , 

UJ is the frequency of the AC field and p is the complex 
dielectric permittivity. Here the subscripts m and p de- 
note the medium and particle. The constants a and b de- 
note half the tube-length and tube-radius. The depolar- 
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On calculating the DEP force of SWNTs using pa- 
rameters from electrical measurements of suspended pure 



SWNtJ*^, the force on m-SWNTs becomes several or- 
ders higher than that on s-SWNTs, indicating that the 
separation should be easy. Furthermore, numerical cal- 
culations suggested that DEP forces of both m-SWNTs 
and s-SWNTs are insensitive to variations in frequency of 
the AC electric field and surfactant concentration within 
the accessible range in experiments^. This does not 
agree with the experimental results where the success of 
DEP separation strongly depended on the frequency and 
the surfactant concentratiorP. Experiments showed that 
the separation was successful only above some certain 
frequencies below which both m-SWNTs and s-SWNTs 
were found on the electrodes. Although the complete 
picture of the transport dynamics is yet to be revealed, 
Krupke et aP attributed the discrepancy to the electrical 
conduction through the surrounding surfactant. They 
have represented the surface conduction effect with the 
effective electrical conductivity which was obtained by 
fitting the solution of Eq. (1) to the experimental results. 

Previous numerical simulations on DEP separation of 
SWNTs have usually considered only the DEP forces 
and the Brownian motion assuming that the bulk flow 
velocity is zero. Although the basic concept of DEP- 
separation in principle is simple, the system involves ef- 
fects which may cause bulk flow motions such as elec- 
troosmosis, thermal convection, electrothermal flow. The 
electroosmosis effect is only important at low frequencies 
(less than 10* Hz, which is well below the frequency ap- 
plied in experiments)*^, and the thermal convection is 
expected to be negligible in micro and nanoscales. On 
the other hand, the electrothermal effect is known to be 
substantial in micro systemS^^EH, Electrothermal flow is 
driven by a body force caused by electric field acting on 
gradients in permittivity and/or conductivity due to a 
non- uniform temperature field*^. 

In this paper, we investigate the impact of the elec- 
trothermal effect on the DEP-separation of SWNTs by 
formulating a dynamical model of the integrated system. 
Here we mainly discuss the transport of s-SWMT, to 
which electrothermal force has non-trivial effects. We 
demonstrate that for a commonly used surfactant and 
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electric field with practical magnitude, electrothermal 
flows can be sufliciently large to have a substantial im- 
pact on the separation efliciency of the DEP method. It is 
shown that electrothermal flows can significantly weaken 
the DEP-separation by driving the s-SWNTs towards the 
electrodes, which is consistent with the experimental ob- 
servations. 



II. MATHEMATICAL MODELING 

Using data for a commonly used material in 
experiment^, the system under consideration consists of 
SWNTs with a diameter of 1.4 nm and a length of 1 
fim dispersed in aqueous solutions of sodium dodecylben- 
zene sulfonate (SDBS). For the electrical conductivity, we 
adopt the effective conductivity (0.35 S/m) suggested by 
Krupke et aP. While the influence of surfactant on the 
electrical conduction of SWNTs surrounded by SDBS is 
not clear, we adopt the empirically effective value. Simi- 
larly, the permittivity of s-SWNTs was set to be 5e^. Al- 
though the magnitude of permittivity of s-SWNTs is ar- 
guable, a variation within one order of magnitude should 
not affect the current analysis. 

Electrical conductivity of SDBS solutions (Tm was es- 
timated as 4 mS/m, 29 mS/m and 230 mS/m for con- 
centration of 0.01%, 0.1% and 1%, respectively^. The 
permittivity of the SDBS solutions em at room tempera- 
ture is about 80eo and has negligible dependence on the 
concentration. It is straightforward to calculate the DEP 
force factor for s-SWNTs and the results show that it de- 
creases about one order when the concentration of SDBS 
increases by one order. 

We solve the heat conduction equation with a heat 
source due to an electric fielcP^, /cV^T -I- a^E'^ — 0, 
where the thermal conductivity fc=0.6 Wm~^K~^. As 
a typical system geometry in experiments^, a pair of 20 
/im-long and 1 /im-wide electrodes were located with a 
distance of 10 /im at the bottom of the calculation do- 
main as shown in Fig. 1(a) The boundary conditions 

dT 

are: Twaii — 300 K at the surrounding walls, and — — — 

dz 

at the bottom, which means the substrate and electrodes 



are thermally insulated. Figure 1(b) shows the tempera- 



ture profiles for various SDBS concentrations above the 
tip of the electrode {x, y)—{0 /zm, 5 (mb), with the applied 
AC potential <j) — 20 V (peak to peak). Maximum tem- 
perature differences are 0.3 K, 1.6 K and 11.8 K on the 
electrode (z=0) for 0.01%, 0.1% and 1% SDBS solutions, 
respectively. 

The electrothermal body force acting on the bulk fluid 
due to the AC electric field is expressed a,^ 



Z (Tn, 



1 



1 + (c^t)2 



EpVem], (2) 



where r — em/cm is the charge relaxation time of the 
solution. The first term on the right hand side expresses 
the Coulomb force and the second the dielectric force. 
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(a)Geometry of electrodes used (b)Temperature changes with z 
in calculations at {x,y)={0 ^m, 5 /jm), </> is 20 

V, a;=300 kHz 



FIG. 1: Geometry and temperature profiles. 



For a typical electrolyte solution, (l/gm) / {dc-ax/ dT) = 
0.02 K-i and {l/e^)/{de^/dT) = -0.004 K"^. The 
first (second) term is dominant in the low (high) fre- 
quency regime (the crossover frequency is uoc = I/t, 
which is 5.6 MHz and 41 MHz for 0.01% and 0.1% SDBS 
concentrations). In the calculation, we adopted 300 kHz 
and 1 GHz, which represent the low and high frequency 
regime, respectively. For low Reynolds numbers, the 
transport term of the Navier-Stokes (N-S) equations can 
be neglected. Hence, the N-S equations become 



9uf 



(3) 



For low Reynolds numbers, the inertial effect can be 
neglected^^. The terminal velocity of s-SWNTs is cal- 
culated as u = Ud + Uf + Ub, here = Fdep/A 
with the translational friction coefficient for a pro- 

67rr;a [r2| 



late eUipsoid given by /t = , /o /^,^ ■ 

/n(2a/o) 

istic variables of electrothermal flow are U 
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L=l /im and t* = 



dT 
dT 



k{l + {LUT)^)r]L' ^ k ' -r^-— 

+ (^''") )'n ^ Here (j) is the applied peak-to-peak AC 
25em0'* 

voltage. From the Einstein relation, the Brownian ve- 
locity can be derived as U{, = ■\/6-D/ dt with the dif- 
fusion coefficient D = {ksT/ft), and dt is the time 
interval of observation. We deflne dimensionless vari- 
ables u = u/U*, f ^ {T - TwaiO/T*, i = x/L, dt = 
dt/ t* and E = EL/cf), dimensionless parameters Pi — 



(l+(^r)2) 
12 



Pa = 0.0463 



dT b^alnmk{l + {u;TY) 



da,. 



and 



^3 = \ a TTT^ 75 • Consequently, 

y ocTm 0.6 • em(p 7ra 

we obtain a set of system equations in a non-dimensional 
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form, 

V-E = 0, V^f =-|E2|, 

^ = V^Uf- Vp+(Vf •E)i^ + Pi|E2|vf', 
ot 

ub = r\/(r + rwaii/ r*)/ 4 Ud = vjEp 

u = Paiid + Uf + P3Ub, (4) 

where f denotes a random number with normal distribu- 
tion whose mean is and standard deviation is 1. Equa- 
tion (4) is solved using finite element method toolbox 
femLegcP^. 

III. FORCE ESTIMATIONS 

Since the inertial effect is neglected, the dimensionless 
velocities are equivalent to the dimensionless forces. Tak- 
ing the characteristic timescale t* as the time interval dt, 
the dependence of the relative magnitudes of DEP, Brow- 
nian and electrothermal forces on the applied electric po- 
tential was examined. Figure [2] shows the dimensionless 
forces at {x,y, z)={0 /zm, 14 /im, 1 ^m) for two different 
surfactant concentrations 0.01% and 0.1%. The y po- 
sition was chosen to be in the regime with the positive 
electrothermal force that attracts SWNTs towards the 
electrode and z=l /im equals to the length of SWNTs. 
Here, the frequency is 300 kHz, which is within the fre- 
quency range commonly used in experiments^. Figure [2] 
shows a clear dependence of the locally dominant force 
on the applied voltage and surfactant concentration. The 
change of the Brownian motion due to voltage and surfac- 
tant concentration is subtle and relatively negligible. The 
voltage-dependence of the strength of the electrothermal 
force is proportional to 0^ and that of the DEP force is 
proportional to (f)'^. As a result, for both concentrations 
0.01% and 0.1%, the electrothermal force overcomes the 
DEP forces at the high voltage limit since the electrother- 
mal force increases with the voltage much faster than the 
DEP force. 

As for the influence of the surfactant concentration, 
the electrothermal force increases with the concentration 
while the DEP force takes the opposite trend. When 
surfactant concentration is 0.01%, the dominant force is 
the Brownian force for low voltage < 6 V), the DEP 
force for intermediate voltage (6 V< < 19 V) and the 
electrothermal force for high voltage {(p > 19 V). On the 
other hand, when the concentration is 0.1%, the local 
DEP force around the selected location is never dominant 
for the entire range of (f>. This implies that for SDBS con- 
centration higher than 0.1%, the positive DEP force on 
s-SWNTs plays minor role and the transport is governed 
mainly by the electrothermal force. The result suggests 
that there is a crossover between the electrothermal and 
DEP force with respect to the SDBS concentration. This 
means that beyond the crossover, an attempt to enhance 
the DEP-separation efficiency by increasing the surfac- 
tant concentration, i.e. by reducing the magnitude of 
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FIG. 2: Forces change with <j) at {x,y,z) — {0 /im, 14 pm, 
1 /Mm). The Brownian force profiles overlap on each other 
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FIG. 3: Forces change with z, {x, y) — {0 /im, 14 /xm), <^=20 V. 
The Brownian force profiles overlap on each other 



the positive DEP on s-SWNTs may result in enhancing 
the electrothermal force, which reduces the efficiency. 

Figure |3] shows the magnitudes of the dimensionless 
forces with the vertical coordinate z for the same (x, y) 
locations as in Fig. |2] Following typical experiments, 
^=20 V was appHed to 0.1% and 0.01% SDBS solutions. 
Figure [3] shows that in 0.1% SDBS, the Brownian mo- 
tion surpasses the DEP force except very close to the 
electrode (less than 1 /im). Here we observe that, in a do- 
main larger than that where the DEP force is dominant, 
s-SWNTs are carried by the electrothermal flow for both 
SDBS concentrations. As a result, the electrothermal 
flow overcomes the Brownian motion in a relatively large 
domain and transport s-SWNTs onto the electrodes. To 
demonstrate this, we will verify the direction of the elec- 
trothermal flow and the actual traces of s-SWNTs under 
the action of the resultant force. 
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IV. 



TRANSPORT SIMULATIONS OF 
SEMICONDUCTING SWNTS 



Figures |4(a) and |4(b) show streamhnes of the elec- 
trothermal flows in 0.1% SDBS solution for two different 
frequencies. They are cross-sectional views of the 3D sim- 
ulation results at the plane of a; = 0. The thick line drawn 
from (x,y)=(0 /xm, 5 /xm) to (0 fim, 25 /zm) marks one 
of the two electrodes [Fig. 1(a) . The evident dependence 



of the flow pattern on cj is highlighted by comparing the 
flow pattern for uj = 300 kHz [Fig. 4(a)| corresponding 
to a typical working frequency with that for lu = I GHz 
[Fig. |4(b)| corresponding to the upper limit of frequency 
explored in experiments. When lu = 300 kHz, the vortex 
core appears close to the gap of the electrodes, whereas 
when w = 1 GHz, it is located close the domain boundary 
{y=25 fim). The pattern of the electrothermal flow varies 
with frequency due to variation in the balance of the 
Coulomb force and dielectric force. The flow patterns for 
u! = 300 kHz and ui = 1 GHz resemble the flow patterns 
for low and high frequency limits, where the Coulomb 
force and dielectric force become dominant, respectively. 
The flow pattern also strongly depends on the geometry 
of the system. In the current 3D system, independently 
of frequency, the electrothermal force gives rise to upward 
flows on the gap of the electrodes and downward flows on 
domain boundary (x = /im, y = 25 /im). This differs 
from the 2D case with an inflnitesimal gap, where the 
direction of the flow circulation at low frequency limit is 
opposite from that at high frequency limiiP^'i^. 

After solving Eq. Q , we obtain the actual trajectories 
of the s-SWNTs as plotted in Fig. [4(c)1 and Fig. [4(d)] The 
simulation times are 33 ms and 24 ms respectively. The 
s-SWNTs are seen to follow the streamlines of the elec- 
trothermal flow, especially well in the region close to the 
electrodes. This agrees with the local force comparisons 
shown in Fig. [3] where the electrothermal force appeared 
to be dominant in a broad region close to the electrode. 
This means that the variation of the flow patterns ob- 
served in Fig. 4(a) and Fig. |4(b)| results in different spa- 
tial trapping distributions of SWNTs. When iv = 300 
kHz, s-SWNTs are attracted broadly on the electrode, 
while when u; — 1 GHz, the s-SWNTs are mainly at- 
tracted to a relatively small region close to the domain 
boundary. 

Let us briefly discuss the influence of the electrother- 
mal effect on m-SWNTs. We have also performed the 
simulations for m-SWNTs with dp = 10^ S/m and — 
— lO'^eJ^. The magnitude of the positive DEP force for 
m-SWNTs is much larger than that for s-SWNTs. Conse- 
quently, the impact of both electrothermal and Brownian 
forces on the overall transport is much smaller for m- 
SWNTs than for s-SWNTs. However, the electrothermal 
force is far from negligible, especially for higher SDBS 
concentrations. In fact, in the gap region of the electrode 
pair where the electrothermal force is directed upwards 
{y=0 /im to 5 fim in Fig. [4]), the electrothermal force over- 
comes the DEP force and repels m-SWNTs away from the 
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(a)Streamlincs of the 
electrothermal flow, 
a;=300 kHz 



(b) Streamlines of the 
electrothermal flow, lli=1 GHz 
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(c)Traces of s-SWNTs, 
4^=300 kHz 
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(d)Traces of s-SWNTs, 
uj=l GHz 



FIG. 4: Streamlines and corresponding traces of s-SWNTs 
in 0.1% SDBS solutions, low (300 kHz) and high (1 GHz) 
frequency regimes with 0=20 V. 



electrodes when the concentration of SDBS is larger than 
0.1%. The impact on the collection of m-SWNTs may be 
limited since the m-SWNTs will be circulated along the 
vortex and eventually transported onto the electrodes af- 
ter a certain time as seen in Fig. |4]for s-SWNTs. How- 
ever, the effect should certainly influence the distribution 
of yielded m-SWNTs on the electrode, which is impor- 
tant for deposition of SWNTs to selected sites for electric 
device constructions. 



V. CONCLUSIONS 

We formulated a dynamical system to simulate the 
DEP motion of s-SWNTs, which is different from pre- 
vious simulations mainly by taking the electrothermal 
flow into account. To realize the DEP separation, the 
applied electric potential needs to be high enough so 
that the DEP force on m-SWNTs overcomes the Brow- 
nian motion. This results in a high temperature gra- 
dient, which creates a substantial electrothermal flow. 
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This electrothermal flow will bring the s-SWNTs close to 
the electrodes where the weaker DEP force on s-SWNTs 
can collect them. Thus, the collaborative action of elec- 
trothermal flow and DEP restricts the allowable range of 
potentials where DEP separation may be possible. The 
higher the concentration of the surfactant is, the stronger 
the electrothermal flow is. Also, when the frequency 
increases, the pattern (direction and magnitude) of the 
electrothermal flow changes. Therefore the electrother- 



mal flow would very likely increase the difliculty both 
in separating s-SWNTs from m-SWNTs and in deposit- 
ing SWNTs to a certain position. The main conclusion 
of this paper is that when designing DEP separation of 
SWNTs, it is necessary to consider the possible influence 
of the electrothermal flow. 
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